Arboreal Bound Entanglement 



Alastair Kay 

Centre for Quantum Computation, DAMTP, Centre for Mathematical Sciences, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, UK and 
Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 

(Dated: August 16, 2010) 

In this paper, we discuss the entanglement properties of graph-diagonal states, with particular 
emphasis on calculating the threshold for the transition between the presence and absence of entan- 
glement (i.e. the separability point). Special consideration is made of the thermal states of trees, 
including the linear cluster state. We characterise the type of entanglement present, and describe 
the optimal entanglement witnesses and their implementation on a quantum computer, up to an ad- 
ditive approximation. In the case of general graphs, we invoke a relation with the partition function 
of the classical Ising model, thereby intimating a connection to computational complexity theoretic 
tasks. Finally, we show that the entanglement is robust to some classes of local perturbations. 



I. INTRODUCTION 

Multipartite entanglement is still a phenomenon that 
is poorly understood and categorised. For some types of 
entangled state, such as iV-qubit GHZ states, it requires 
very little noise (loss of a single qubit) to entirely destroy 
the entanglement, whereas others are much more robust, 
such as the multipartite states based on error correcting 
codes. One is therefore prompted to ask about the types 
of multi-particle entanglement that arise in nature, and 
query how persistent entanglement is within the thermal 
states of local Hamiltonians. The ability of entanglement 
to persist at high temperatures will be crucial to future 
experiments and could have a direct bearing on the con- 
struction of quantum memories [T] . 

We defer the analysis of the persistence of entangle- 
ment in general Hamiltonians to future studies. In this 
paper, we instead concentrate on a specific class, graph 
Hamiltonians, with the intent of gaining insights for these 
future analyses. Graph states are particularly interest- 
ing because they frequently arise within the study of 
quantum information. Many states such as Bell states, 
GHZ states and CSS error correcting codes are all equiva- 
lent, under the action of local unitary rotations, to graph 
states. It is already known that all thermal graph states 
can be distilled up to a finite temperature [5] and that for 
some subset, which includes many interesting examples 
such as the graphs corresponding to square lattices of 
arbitrary dimension, this temperature is tight i.e. above 
that temperature, entanglement can't be distilled [S] [3J . 
On the other hand, for some models, it has already been 
shown that there is still entanglement present in these 
models. This was first demonstrated for GHZ Hamilto- 
nians [3], but has since been shown to exist in simpler 
models, such as for square lattices [5]. 

Papers such as [5j simply give an existence proof for 
bound entanglement, which comes as little surprise since 
it would be quite remarkable if the distillability tempera- 
ture (upper bounded by the temperature at which there 
exists a bipartition of the qubits such that the state is 
positive under the partial transpose operation) coincided 
with the full separability temperature (lower bounded by 



the temperature at which there exists a bipartition of the 
qubits such that the state is negative under the partial 
transpose operation), and do not tackle the question the 
maximum temperature to which the entanglement per- 
sists. This is part of what we address here, although 
we also discuss other issues such as those arising in [5] 
about when one can give a fully separable decomposition 
of the thermal state, as well as giving the optimal entan- 
glement witnesses (in Sec. IV). In Sec. Ill we will analyse 
the linear graph (also known as the ID cluster state) in 
particular detail, and prove that the state becomes fully 
separable as soon as all possible bipartitions are positive 
under the partial transpose operation (PPT), as is the 
case for the thermal states of all tree graphs. In conjunc- 
tion with the previous distillability result, we thus know 
the entire entanglement structure of these states. We 
will demonstrate what features of the proof are particu- 
lar to thermal states by contrasting with the global and 
local depolarising noise cases in the appendices, which 
also serves to resolve some outstanding questions from 
[3J. The techniques that appear here are a generalisa- 
tion of those presented in [7], in which the geometry of 
the star graph was found to be particularly beneficial for 
proving many of the same properties for GHZ-diagonal 
states. For treatment of that special case, we refer the 
reader back to that paper; we will not dwell on it here. 

The treatment of entanglement witnesses in this paper 
is important in its own right. There have been numer- 
ous studies of the entanglement in graph states [8HT3]. 
Many of these detect entanglement on an ad hoc basis, 
i.e. an operator is found which detects some entangle- 
ment. Our strategy rather takes the opposite approach; 
we start from determining the absolute limits of the pres- 
ence of entanglement, describe the witnesses that work 
at these limits, and then show how to approximate these 
on a quantum computer. These witnesses can therefore 
be expected to be substantially stronger than existing 
witnesses in certain regimes (particularly in the case of 
thermal noise), although in some instances these other 
witnesses can also be shown to be optimal [8] [14] . Other 
recent studies [HJ [16] provide an alternative perspective 
- given a set of measurement data, they solve a semi- 
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definite program in order to determine the most and least 
entangled states compatible with those measurement re- 
sults with respect to a specific entanglement measure. 

A. Multipartite Entanglement Structure 

Beyond a simple statement of whether a particular 
graph-diagonal state is entangled, we wish to be able to 
categorise this entanglement. A first natural question is 
how useful is this entanglement? The most useful type 
would be distillable entanglement, and the distillability 
regime has already been investigated 2,3. Beyond that, 
[17j suggested a hierarchy of separability criteria: for a 
state of N qubits, with respect to how many independent 
groupings of the qubits, k, is the state separable? It is 
said to be fc-separable. For instance, if a state can be 
written as 

E1,2...N A _ N A + 1,...N B ^ Nb + 1-..-Nc 
PtPi ^ Pi ®Pi 

i 

for some grouping of the qubits, the state is 3-separable. 
If a state is not /c-separable with respect to some parti- 
tion, we say that there is some fc-partite bound entan- 
glement present (the smaller k, the stronger one might 
consider the entanglement to be). The particular case 
which will arise frequently is when the state is not distil- 
lable, and yet not every bipartition is separable, so there 
exists some bipartite bound entanglement. We empha- 
sise that this is not the same as the genuine multipartite 
entanglement that is often discussed, particularly in the 
pure state case. For instance, in that case, a 3-qubit 
state (| 00) + 1 11)) |0) /V2 is bipartite entangled, whereas 
(|000) + |lll))/v / 2 is tripartite entangled. Since this class 
of entanglement will never arise in the present treatment, 
we are not unduly worried about possible confusion. 

B. Graph State Basics 

Consider a graph G which is composed of edges E and 
N vertices V. With each vertex we can associate a qubit. 
The stabilizer operators are defined as 

K n = X n J] Z m VneV 

{n,m}£E 

and the Hamiltonian for the graph G is 

where X n is the Pauli X matrix applied to the qubit 
on vertex n. We use Z x to denote the application of Z 
rotations to all qubits for which the iV-bit string £ is 1, 
i.e. if we denote the n th bit oi x £ {0, 1}^ as x n , we have 

z*= n z n- 

nev 



K x is similarly defined in terms of the elementary K n . 
The eigenstates of the Hamiltonian are \ip x ) = Z x \ip), 
where denotes the ground state of the Hamiltonian 
(known as 'the' graph state for the graph G). These 
states are all simultaneous eigenstates of the stabilizers 

K n \1> a ) = (-If" 

and they have energies —NA/2 + Aw x , where w x is the 
Hamming weight of x. A general graph-diagonal state p 
(meaning that (tp x \ P \^y) = if x 7^ y) can be written as 

P = 2* ^ SyKy ^ 
ye{o,i} N 

for real coefficients s y with sq = 1 and 

(-!)"% >° Vxe{0,l} N . 

ye{o,i} N 

Any state can be made graph state diagonal via stochas- 
tic local operations while maintaining the values of s y 
[HI [19]. The thermal state is written more succinctly 
since s y = s Wy where s — tanh(/3A/2) and j3 is the stan- 
dard inverse temperature, 

P = TrU-PH) = W II C 1 + tanh(/?A/2)X„). 

This state is the ground state with Z rotations ap- 
plied independently at each site with probability 

_ e-? A 
P ~ 1 + e-^ ' 

This interpretation shows that there is a single critical 
temperature for the persistence of entanglement - it is 
not the case that as the temperature is varied the system 
starts entangled, becomes separable, and then becomes 
entangled again because, through local operations (prob- 
abilistic application of Z rotations), we can convert any 
thermal graph state into one of a higher temperature, 
and this process must only reduce the entanglement. It 
also makes this model of interest in a variety of experi- 
mental scenarios, since the thermal state corresponds to 
local dephasing noise acting on \ip). 

Alternatively, \ip) can be defined constructively. It 
is formed by preparing each qubit in the state |+) 
and applying a controlled-phase gate between every 
nearest-neighbour pair (i.e. every two vertices forming 
an edge). Perhaps the most vital observation, which al- 
lowed the previous proof of distillability [H [3] , is that the 
controlled-phase gate commutes with the Z noise. This 
allows us to map between different graph states without 
propagating the noise in a bad way. 

One can apply local unitary operations to graph states 
and this can map them into other interesting states, such 
as those related to error correcting codes [2Q1E1] . Certain 
local unitaries, however, can map graph states into dif- 
ferent graph states [22] . Such operations will keep graph 
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diagonal states as graph diagonal states (but for a dif- 
ferent graph). However, the permutation of the diagonal 
elements means that these local rotations do not map 
thermal states into thermal states. 

It will often be the case that we are interested in 
two-colourable graphs. In the literature on graph the- 
ory, the term 'bipartite' is more commonly used. How- 
ever, we avoid this terminology in order to avoid confu- 
sion since bipartite will arise in other contexts. For the 
graph-theoretic terminology not defined here, we refer 
the reader to |2"3"1. 



Proof. Z-measurements (which serve to remove vertices 
from a graph state) commute with the local Z noise 
present in the thermal state, so by applying these mea- 
surements (and applying compensating Z rotations on 
the neighbours of any spins we measure in the |1) state), 
we can cut out any spins we want from G. In partic- 
ular, we can remove all the vertices in G but not in H. 
This leaves behind H , from which we conclude that if en- 
tanglement persists in any bipartition of H , it must also 
persist in an identical partition of G (where it is irrele- 
vant on which side of the partition we place the additional 
spins) . □ 



II. BIPARTITE BOUND ENTANGLEMENT 

In [21 [3] , we described a class of graphs whose thermal 
states can be optimally distilled. Every graph state can 
be distilled up to the temperature 



T, r , 



A 



distillablc 



k B ln(V2 + l) 



Clearly, it is necessary for distillability that there should 
be distillable entanglement between every possible bipar- 
tition, i.e. every bipartition should be non-positive with 
respect to the partial transpose operation (NPT). 

Our program now consists of deriving when there ex- 
ists at least one NPT bipartition of a thermal graph state, 
and comparing it with the temperature for distillability. 
Let us emphasise that there could still be entanglement 
above this bipartite bound, which would correspond ei- 
ther to bipartite entanglement not detected by the NPT 
condition, or to genuine multipartite bound entangle- 
ment, but we need to establish this bound first. We start 
from Ecm. ([I]) and introduce a bipartition z g {0, 1}^ 
i.e. all the vertices with z n = are on one side of the 
partition and those with z n — 1 are on the other side. 
We will take the partial transpose on the z n = 1 side. 
Recall that under the partial transpose, the Pauli oper- 
ators alter by Z n M- Z n , X n i— > X n but Y n n- (— l) z ™Y n . 
Thus, 



P 



1 



E 

y£{o,i} N 



Observe that products of stabilizers remain as products 
of stabilizers and, as a result, the eigenvectors of p PT are 
just \tp x ), with eigenvalues 



1 



E (-!)*'%(- 

ye{Q,i} N 



For thermal states, we are interested in finding the small- 
est positive value of s such that the smallest eigenvalue 
is zero. In the remainder of this section, we will establish 
the optimal choices of both x and z. 

Lemma 1. Let G be a graph with an induced subgraph, 
H . Entanglement persists in the thermal graph state of 
graph G up to at least the level it persists for H . 



We learn from Lemma [T] that bound entanglement is 
more persistent in large systems than in small ones! 

Let us neglect the factor of 2 , since it is of no conse- 
quence. The eigenvalues are thus denoted by 



/£(«)= E (-ir^(-i)^.w^^ 

y&{o,i} N 



») 



for the thermal state. By neglecting the subscript x, we 
mean the case of x = {11. ..1}, and by neglecting z 
we mean (for two-colourable graphs) the two-colouring 
partition. 



f N (s)= E 



(2) 



We will be particularly loose with the superscript N. As 
an N, it denotes the number of vertices in a graph (and 
hence the number of bits of x and z) , which is useful for 
recursion relations. However, we might instead use G to 
denote a particular graph and, often, it can be neglected 
entirely as it is clear from the context. 

The notation of neglecting the x and z subscripts will 
arise because we can prove the best choices: 

Lemma 2. For thermal graph states, the smallest pos- 
itive value of s for which min X)Z /-^(s) = is achieved 
with x — {11 ... 1}. 

Proof. Consider a bipartition z. Select a spin n, using 
x to index all bit strings on the N — 1 bits not n, and 
use k as an N — 1 bit string which is except in the 
positions which have an edge to spin n which crosses the 
bipartition, k m = E m , n {z m © z n ). One can show that 



/x G o (*)-/£(*) 



(3) 



From Lemma [T] we know that the graph G is at least as 
entangled about any bipartition as it is when qubit n is 
removed. Therefore, in the region when all the f^ z (s) 
are positive (i.e. for values of s > but smaller than the 
critical value), all the /^"(s) are also positive. There- 
fore, in this region, the smallest of all the f& z (s) is for 
x = {11 ... 1}, i.e. this is the one that will become 
first. □ 
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Lemma 3. The optimal bipartition for the thermal state 
of a two-colourable graph is a partition into its two colour 
classes. 

Proof. Consider each pair of nearest-neighbour variables 
z n © z m as if it were an independent variable z nm . We 
can now apply the same technique, comparing when one 
of these variables is either or 1: 



/ x,(z\nm)0 ( S )~fx,(z\nm)l( S ) 



2 f G\r 



(x\n,m)(&k,z\(nl) , (ml) 



where the bits of k, ki — E ni iz n i + E mt iz m l- Again, in the 
region we're interested in, the right hand side is positive, 
meaning that the optimal choice is to set the value z nm = 
1. Thus, we should set all pairs z n © z m — 1. In general 
this is not possible because the variables are not really 
independent. However, for two-colourable graphs, the 
natural bipartition allows us to do this. □ 



A. The Classical Ising Model and Computational 
Complexity Considerations 

We want to evaluate f x ,z(s) for some bipartition z of a 
thermal graph. With a fixed bipartition, then the partial 
transpose condition is unaffected by unitaries that are 
local with respect to the bipartition. As such, we can use 
controllcd-phase gates to remove any edges of the graph 
that are on one side of the bipartition, and this always 
reduces us to a two-colourable graph on which we wish 
to evaluate Eqn. ^ and for which z is the two-colouring 
bipartition. We could just as easily rewrite this as 



/(*) 



E 

ye{o,i} N 



which is equivalent to the partition function of the clas- 
sical Ising model for the same graph, 



P '{T, n h n S n +J E { „, m}e B S„S m +kl) 



z = E ' 

se{-i,i} 



where, since the S n take values ±1, we rescale the y n with 
S n = 2y n — 1. To complete this mapping, we identify 

0J = in/A 

(3'k = K H - a ) + ™\E\ 

p'h n = §ln(-s) + ~d n , 

where \E\ corresponds to the number of edges in G and 
d n is the degree (coordination number) of site n. Of 
particular interest would be a solution in 2D. However, 
the general solution to the 2D Ising model for arbitrary 
h is unknown. Indeed, there are known instances of field 
strengths for which the problem is NP-hard [24 . This 
leads to the suspicion that identifying the exact critical 
temperatures is a hard problem. 



From [33] we know that /(l) can be efficiently 
evaluated. However, were we able to evaluate f(s) at 
several values of s, we would be able to resolve the values 
of 

(_iyy T ltr(A)y 

ye{o,i} JV 

Wy—k 

for all k, where ltr(j4) is the lower triangular compo- 
nent of the adjacency matrix A (this is equivalent to the 
graph theoretic problem of calculating the number of in- 
duced subgraphs of G of k vertices which have an odd 
number of edges) . There is a closely associated problem 
for which there exists a dichotomy theorem [37] between 
instances of the problem which are efficiently solvable 
and those which are #P-complete, based on whether the 
logical expressions represented by . ltr(A) i jy i yj can 
be written using a constraint language which is affine of 
width 2 (i.e. as a string of conjunctions, where each term 
is only an XOR of two variables g/j), and there are many 
#P-complete tasks that remain so, even under the re- 
striction to planar bipartite graphs [28j . Unfortunately, 
however, classifying the graphs that would give us an effi- 
cient solution is related to a long standing open problem 
in complexity theory, ©2sat [55]. As such, while we are 
left believing that evaluating f(s) is hard at most values 
of s for certain graphs, we have not succeeded in proving 
this. 



III. ID CLUSTER STATE 

There are many graphs that we could consider, and try 
to evaluate the critical temperatures. Since our study is 
motivated by the wish to understand the entanglement 
properties in thermal states of naturally arising Hamilto- 
nians, we ideally wish to understand regular lattices such 
as the ID and 2D square lattices, and compare them to 
the case of GHZ states studied in [7J: does the finite de- 
gree of each vertex mean that entanglement only persists 
to a finite temperature? We shall start by applying our 
formalism to the ID chain. We want to determine when 
f N (s) = in order to find if the state is entangled ac- 
cording to the PPT criterion. If f(s) is negative, the 
state is certainly entangled, but, if the value is positive, 
we may learn very little. One can verify that 



/" 00 = (1 + -2«/"-'( fl ), 

which can be solved to find f N (s) for a specific s. The 
solution is of the form f N = Ar+ + Br^ , where r± are 
the roots of the equation 

r 2 - (I + s)r + 2s = 0. 

The coefficients A and B allow us to match the initial 
conditions of / 1 (s) = 1 - s and f 2 (s) = 1 - 2s - s 2 . 
We want to find the value of s for which f N (s) —> in 
the large N limit. This is readily achieved if |r| < 1. 
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FIG. 1: The critical temperature for the existence of bound 
entanglement in the linear cluster state compared to the re- 
gion of distillability and the N — > oo limit. 



However, we also want to know the largest value of s 
that achieves this without undergoing a sign change as N 
increases - this assures us that we have the smallest root. 
These sign changes are brought about only if the r± are 
complex. The roots are real in the region < s < 3—2^/2 
and thus we have f N ^(3 - 2^2) -> without any sign 
changes, which means that in the thermodynamic limit, 
NPT bound entanglement persists until tanh(/3A/2) = 
3 — 2\/2. Fig. [I] compares this limit to the numerical 
values calculated for finite N. 



A. Full Separability 

We now know a temperature up to which there is cer- 
tainly a bipartition with respect to which the system 
remains NPT and hence entangled. It could be that 
NPT doesn't detect all the bipartite bound entanglement 
or that there is some multipartite bound entanglement 
present. Hence, it will be useful to establish a limit on 
when there can be entanglement by giving a threshold 
temperature above which the state is fully separable. If 
this comes out to be the same as the above PPT limit, all 
the entanglement is bipartite bound entanglement, and 
is detected by PPT. 

We use essentially the same considerations as in [7], 
that is to say that we define two stabilizers K x and K v to 
be compatible if they have a simultaneous product state 
decomposition. A sufficient condition for this is that for 
any site where K x is a Pauli matrix a ^ 1, at the same 
site, K y should either by a or 1, and vice versa 1 . For 
instance, on a 3-qubit chain, the 3 terms 

K 1 = X®Z®1 
K 3 = 1 ® Z ® X 
K t K 3 = X (g) I® X 



This is not a necessary condition and, for completeness, we will 
illustrate an alternative way of making the decomposition in Ap- 
pendix [D] 



have a common product basis of |±) |0) |±) and 
|±) |1) |±), but are incompatible with K%. For a 3-qubit 
chain, these are the only compatible terms. This means 
that the expression 

1 + ^1+^3 + K t K a 

is a separable state (we have added sufficient 1 such that 
the smallest eigenvalue is 0). So, let us return to the 
3-qubit thermal state 

f 3 

n=l 

and expand it as 

8p = s 2 (l + Ki)(l+K 3 ) 

+s{t + K 2 ) + s 2 {t + K 1 K 2 ) 
+s 2 (l + K 2 K 3 ) + s 3 (l + KxK 2 K 3 ) 
+(s ~ s 2 )(l + K\) + (s - s 2 )(l + K 3 ) 
+1(1 - 2(s - s 2 ) - s - 3s 2 - s 3 ). 

By construction, each of the terms is a separable den- 
sity matrix provided the coefficients in front of them are 
positive. In particular, this means that the coefficient in 
front of the 1 term, < 1 - 3s - s 2 - s 3 = f 3 (s) . 

We now aim to perform a similar construction for 
chains of length N. We will define a basic string y to 
be a string which is of the form 00 . . . Olf . . . f 00 ... 0, i.e. 
it is a continuous block of w y Is, with some unspecified 
number of 0s on either side. A basic stabilizer is a stabi- 
lizer K y where y is a basic string. Two basic stabilizers 
have a compatible basis if and only if their basic strings 
are non-overlapping, and separated by at least one ver- 
tex. Let by denote the decomposition of y into its set 
of compatible basic strings, i.e. K y = Y\ x£b K x . The 
members of b y are the generators of a group B y . We 
are also interested in the set of strings C y defined such 
that x € C y if y € B x , i.e. it is the set of strings which 
are compatible with y and have all bits set to f if the 
corresponding bit of y is set to 1. Since a basic string y 
necessarily satisfies the property 

(_!)«» (-i)!^ 1 »»»»+! = _i, 

then evaluation of ( — l) Wy ( — f )2„=i for any gen- 

eral y simply returns the parity of \b y \, the number of 
basic strings that y is formed from. 

With these definitions in place, it is easier to state what 
we wish to prove, 

V x£b y 

where 

f y (s) = (-l)IM J2 s IUj; (-1) IM - 

X£Cy 
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This quantity f v {s) is exactly the same as evaluating 
s Wy f{s) over the induced subgraph specified by taking 
the chain and removing all the vertices, and their neigh- 
bours, where y n — 1. Hence, we learnt from Lemma [T] 
that in the region where the state is PPT with respect 
to all bipartitions, all f y {s) > 0, and the interpretation 
of the decomposition giving a separable state will be cor- 
rect. 

In order to prove our desired decomposition, we start 
by reverting to the more general notation of 

P = 2W E] s yKy, 
v 

which will make it easier to keep track of a given term 
s y = s Wy . Observe that if C y = y, i.e. y is not compati- 
ble with any strings of weight greater than w y , then the 
coefficient f y (s) = s y . 

If y is not a basic string, K y will decompose, and we 
have to remove s y from the coefficients in front of all the 
K x terms where x G B y . If K y = K x K a and a is a ba- 
sic string, then the coefficient in front of the K x term 
will be altered by —s y because of the weight s y K x in the 
s y Y\ zeby {^ + Kg) expansion. Now consider instead that 
K y = K x K a Kb where a and b are basic strings. The 
expansion of K y has altered the K a K x and K^K X coeffi- 
cients by — s y , and the expansion of all 3 of these terms 
alters the K x coefficient by 2s y — s y = s y (in addition 
to the s a+x and Sb +X coefficients which we can follow 
through independently in exactly the same way). It now 
follows inductively that if K x and K y differ by r basic 
strings, the K x coefficient is altered by (— l) r s y . Now, 
from our previous observations, we have that 

(— l) r = (_i)IM+IM. 

So, if we start with p = Yly s yK y working from the 
strings y of highest weight, replacing K y with fl^efe,, (1 + 
K x ), we end up with the decomposition in Eqn. Q, 
which is separable provided all the coefficients are posi- 
tive. 

It follows immediately that for all thermal graph states 
where the graph is a chain, or indeed a tree, the state be- 
ing PPT with respect to all possible bipartitions is nec- 
essary and sufficient for separability. 

Why does this proof apply to trees but not to more 
general graphs? It was vitally important in the con- 
struction that the compatibility of the stabilizer oper- 
ators was hierarchical, meaning that for every x € B yi 
B x C B y , In a ring of N qubits, for instance, this does 

not hold - Yln=i K n is compatible with n«=i Km and 
Hn=0 K ^n+i, an d then Ilnio^ri+i is compatible with 
each iC2n+ij but rXn=i is n °t compatible with K n . 
In Appendix [C] we will show how far we have been able 
to extend our results to this more general case. Can the 
proof extend to graph-diagonal states other than ther- 
mal states? There will certainly be some for which this 
is true, and in Appendices [X] and [B] we look at some 
common noise models. 



IV. ENTANGLEMENT WITNESSES 
SATURATING PPT 

Using the formalism that we've developed, it's rela- 
tively simple to find an entanglement witness for a given 
graph G that will saturate the PPT threshold for any 
state which is diagonal in the graph state basis. To do 
this, we measure the observables 

w *,*=^w E (-ir y (-i)^^ ynvm(zn(Bzm) Ky, 

ye{o,i} N 

For any arbitrary density matrix p, we calculate 
ye{o,i}« 

i-e. f x , z /2 N , the eigenvalues of the partial transpose of 
the state about bipartition z. Hence, for a graph diago- 
nal state, finding Tr(Wx, z p) < for any x or z proves it's 
entangled. Moreover, such witnesses saturate the PPT 
bounds for this class of states. An important observa- 
tion is that this is a genuine entanglement witness i.e. for 
any state p — J2 X y Px,y |A X ) (Aj,|, which may not be di- 
agonal in the graph state basis, finding one of the observ- 
able to be negative witnesses the fact that it's entangled. 
To prove this, note that any p can be converted, via lo- 
cal probabilistic operations, into a graph diagonal state 
Pd = J2 x Px,x \^x) (A^l with the same diagonal elements 
[15] . and hence the same values of s y — J2 X Px,x{— 
So, if p is fully separable, it will have the same value 
of Ti(Wp) as pd, which we know will be positive since 
the local conversion to a diagonal state cannot introduce 
entanglement. 

Since trees have the existence of an NPT bipartition 
as a necessary and sufficient condition for the existence 
of entanglement in the thermal state, the witnesses are 
optimal in this case. 



A. Implementing the Entanglement Witnesses 

So far in this section, we have seen how entanglement 
witnesses can be designed for any given graph which work 
up to the transition between NPT and PPT. For a vast 
range of GHZ states [7], and the thermal states of trees, 
we know that this transition point is also the transition 
to full separability. Thus, in experiments, we will be ex- 
tremely interested in measuring these witnesses in order 
to detect entanglement. However, measuring Tr(W^; )Z p) 
appears to be a daunting task since W XfZ is a sum of 
exponentially many terms K y . We will now show how 
this worry can be circumvented; we can efficiently mea- 
sure Tr(W XjZ p) up to an additive approximation. The 
protocol that we will follow is largely inspired by studies 
of the evaluation of tensor networks on a quantum com- 
puter |30j which proved particularly useful for calculating 



7 



partition functions in classical statistical mechanics. Re- 
lated problems also arose in |31j . 

Let us assume that we are given k copies of p with 
which we are to estimate Tr(W XlZ p). This is equivalent 
to combining the probabilities of detection of each of the 
graph states \ip a ) via 

X a(lpa\p\lpa) , 

a£{0,l} N 



where 



2 N 



The na'ive approach would be to measure each of these 
probabilities independently, and then combine them, but 
that is a very inefficient approach. Instead, we would like 
to measure Tr(Wp) directly. Consider the circuit 



\0)a 
P — 

1/2 N 



H 



CP F _ 



H 



Z 



H 



CPi 



The top wire, A, is a single ancilla qubit, whereas the 
other two wires represent N qubits on which gates are ap- 
plied transversally, except for CPe which indicates that 
controlled-phases need to be applied between the qubits 
of a given register for those pairs of qubits corresponding 
to an edge of the graph. CP E is a modification of this 
in which a controlled phase is only applied between two 
vertices n and m if joined by an edge of the graph and 

For any p — ^ a b pt.a | A&) (A a |, the first two gates are 
just the sequence that maps |A a ) into a computational 
basis state |a), so this circuit simply represents the input 
of a state 



° = /^Pb,a \b) (a\ 



a.b 



i 



to a Hadamard test using a 2N qubit unitary V. The 
Hadamard test has a probability of finding the ancilla in 
state |0) of 

Po = ±(l+Tr(<xF)) 

if Tr(aV) is real. After k repetitions, the probability of 
incorrectly estimating the value of Tr(aV) to an accu- 
racy e is bounded by Chernoff's bound to be no worse 
than 2e~ 2ke , so using k ~ 0(e~ 2 ) gives a constant fail- 
ure probability, independent of N. It remains to argue 
that Tt(ctV) = Tr(W x ,zp), which is readily verified. We 
conclude that using k <~ 0(e~ 2 ) copies of p allows us to 
estimate Tr(Wp) (which is bounded between ±1) to an 
accuracy e, independent of system size, with a number 



of gates O(Nd) where d is the maximal degree of any 
vertex. 

The difficulty with any additive approximation is the 
accuracy scale that we need to achieve. Given that the 
eigenvalues are bounded between ±1, we don't expect 
to significantly outperform the protocol we've just de- 
scribed. However, if we look at the ID thermal chain, for 
instance, then for the two-colouring bipartition and x 



{11 ... 1}, the size of the eigenvalue scales as — (is 



,N/2 

(see Appendix |AJ , which shows that we would require 
0({2/s) N ) copies of p to have any hope to detect entan- 
glement using this bipartition, and we will most likely 
be doing this at finite s. How, then, are we to proceed? 
There is one saving grace in that we can sum many eigen- 
values together. Indeed, this happens naturally in our 
protocol as a result of the trace operation. Clearly, if 
several eigenvalues sum to give a negative value, then at 
least one of them is negative, and so this also acts as an 
entanglement witness. For instance, we can use the bi- 
partition z to define a subset of vertices, S, which have 
at least one edge crossing the bipartition. The vertices in 
S' = V \ S can then be disentangled from the graph state 
without changing the eigenvalues of the partial transpose 
operation. If we just trace out those qubits, however 
(without applying the transversal controlled-controlled- 
phase on those qubits), then we end up summing 2' s > 
eigenvalues. For thermal states, this provides the abil- 
ity to measure the entanglement witness on any induced 
subgraph of \S\ qubits, and the required accuracy is only 
exponential in | *S" | not N. The reason for it being partic- 
ularly useful for the thermal state is that all of the qubits 
that we remove are separable from the rest of the system, 
in the state |(l + s) |0) (0| + A(l-s) |1) (1|. Since for the 
thermal state (see Lemma [T]) these bipartitions/induced 
subgraphs induce a partial ordering in the entanglement 
detection ability, each of which can detect entanglement 
up to some finite temperature, this is extremely promis- 
ing for the detection of entanglement in this case. More- 
over, we observe from Fig.[T]that even short lengths (15*1 
independent of N) of induced subgraph can detect almost 
all of the temperature range of any longer chain. 



1. Evaluating Critical Temperatures on a Quantum 
Computer 

The technique that we have used for implementing the 
entanglement witnesses also suggests that we might be 
able to implement a quantum computation to evaluate 
f(s) to some level of approximation, such that we might 
search for the critical temperature. Although the prob- 
lem of exponentially vanishing eigenvalues will remain, 
we will now see how the exponent is at least reduced 
compared to the witnesses above. 

Consider the following circuit 



10), 



H 



H 



where 



3.0 



P = 



+ 



(|0) (0|+a|l)<l|) { 



The use of p is what distinguishes this algorithm from 
the entanglement witnesses - if we are provided with the 
thermal state, there is no simple way to convert it into 
this form of p, but we can certainly create it if we know 
the value of s that we want to test. The probability of 
getting the |0) result on the top register is given by 



1 



/(*) 



(1 + *) 



N 



which allows us to evaluate f(s)/2 N (the eigenvalues of 
the partial transpose operation), although exponential 
accuracy is still required in order to detect entanglement 
close to the separability threshold. There are some fur- 
ther improvements that can be made using the obser- 
vations made in App. |c] for the classical algorithm 2 , al- 
though they do not significantly impact the scaling. 



V. STABILITY TO PERTURBATIONS 

Calculating the temperature at which a particular 
thermal state becomes PPT with respect to all bipar- 
titions, and being able to detect it, are important steps. 
In the real world, however, we will never generate exactly 
the Hamiltonian that we want. Thus, we must examine 
how stable the entanglement is to small perturbations. 
The simplest possible case is to redefine the Hamiltonian 
as 



H 



I5> 



with s n = tanh(/3A„/2), corresponding to different de- 
phasing rates on each qubit. Due to the numerical work 
of [S] , we already expect that the critical temperature will 
not be significantly affected. In the case of a ID chain, 
we can derive a recursion relation for the polynomial, 

f(s k . . . s N ) = (l + s fe )/(s fe+1 . . . s N )~2s k f(s k+2 . . . s N ), 

which yields the temperature at which the state becomes 
positive with respect to the partial transpose taken on 
any possible bipartition. Numerically, we could sample 
the A„ from some distribution - we chose a Gaussian 
distribution of width a - and calculate the critical tem- 
perature. This has been plotted in Fig. [2] indicating a 
substantial robustness. Analytically, we can note that 
mm Sk = 3 - 2V2 gives a lower bound on the critical 
temperature in the thermodynamic limit. 



2 Namely, one can use specifics of the graph structure, such as 
its bipartite nature, and the existence of a subgroup within the 
group of bit strings {0, 1}^ under addition modulo 2. 
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FIG. 2: Comparison of the distribution of critical tempera- 
tures (the error bars) with the unperturbed critical temper- 
ature (the circles) for a thermal chain when the weights of 
the stabilizers are perturbed according to a Gaussian distri- 
bution, centred on A = 1 with a = 0.1, and then averaged 
over 100 samples. 

While this is a relevant perturbation to consider, per- 
turbations are likely to exhibit much greater variety. 
There is a further case where we can make significant 
progress, although it remains far from the case of arbi- 
trary local magnetic field: we add a perturbative field of 
the form 

JV 

-I V A 7 

n=l 

As \5\ increases, the ground state is progressively less 
entangled. The terms K' n = A n K n + S n Z n mutually 
commute, [K' n ,K' m ] = 0, and K' 2 = (A 2 + 5 2 )t, so one 
can write the thermal state as 

, N 

P = W n { l + tanh^TAlT^iC/VA 2 + <P 

n=l 

For simplicity of notation, we'll take all the A„ equal, and 
all the 5 n equal, although none of the subsequent analysis 
is contingent upon it. The thermal state is expanded as 



E E 

x£{0,l} N yG{0,l} N 



Va 2 + s 2 



Z-rK,. 



although the summation over y is restricted to cases 
where y n = if x n = 1. We can now take a partial 
transpose with respect to a bipartition z. 



2 N p PT 



E E 

^£{0,1}™ ye{o,i} N 



n/A 2 



S 2 , 



No longer is this a product of commuting terms, so the 
eigenvectors are not as obvious as the case where 5 = 0. 
Nevertheless, the minimum eigenvalue is no larger than 
(ip\ p PT \(p) for any \ip). If we take \tp) = \ip x ), this is 
entirely equivalent to measuring the expectation value 
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of the entanglement witness W x>z . Let us take \ipn...i) 
This gives 



E 



2/6{0,l}« 



-sA 



VA 2 + S 2 



since 



Ti{K y K z Z x ) = 2 N 5 y , z 5 i 



y,z<J x flO...O- 

This is exactly the expression we had for no perturbation, 
except that the effective s has changed. The critical tem- 
perature 0$ at which the expectation value of this state 
is zero is hence related to the unperturbed /?o by 

A 



VA 2 + S 2 



tanh(i/WA 2 + <5 2 ) = tanh(§A,A). 



Furthermore, j3$ is an upper bound on the true critical f3. 
If one examines this expression, it is remarkably robust. 
For the infinite 1-D chain, we find that the lower bound 
on the critical temperature T decreases monotonically, 
but is still at 99.0% of the unperturbed value for 6 = 1. 

In the special case of a ID geometry, we could make 
progress in analysing other perturbative fields such as X, 
Y and ZZ through transformation into a standard form 
which can be analysed via the Jordan- Wigncr transfor- 
mation [32 [33] . However, we have not succeeded in an- 
alytically calculating the expectations W XyZ . 



VI. CONCLUSIONS 

We have shown that, for the ID thermal graph state, 
bound entanglement persists to a finite temperature, in 
contrast to the extremely non-local GHZ systems [7], 
which have entanglement persistence that grows with sys- 
tem size. In the thermodynamic limit, the critical tem- 
perature for the chain is 



A 



separable 



k B InvT 



in comparison to the distillability threshold 



T, 



A 



distillablc 



k B ln(V2+l) 



Note that as the length of the chain increases, so does the 
temperature for the persistence of bound entanglement, 
unlike the distillable entanglement. Furthermore, this 
entanglement is all bipartite bound entanglement; there 
is no multipartite bound entanglement present. 

We have given entanglement witnesses that detect all 
the PPT entanglement that is present in any graph- 
diagonal state of any graph, and are hence optimal if 
the state is the thermal state of a tree - the state is fully 
separable if not NPT. We have discussed the implemen- 
tation of these witnesses on a quantum computer. The 



entanglement in the thermal state is robust to pertur- 
bations in the Hamiltonian. In the appendices, we have 
discussed how optimality of the PPT condition might 
extend to other graph-diagonal states for chains, and for 
thermal states of graphs containing loops. We conjecture 
that for fc-colourable graphs, the thermal state becomes 
fully separable when it becomes fc-separable. However, it 
seems likely that determining the critical temperature for 
graphs such as the 2D lattice will be a computationally 
hard problem. Resolving these questions are interesting 
avenues for future study. 

This work was supported by the National Research 
Foundation & Ministry of Education, Singapore, and 
Clare College, Cambridge. The author thanks Simone 
Severini and Ravishankar Ramanathan for useful con- 
versations. 



Appendix A: Global Depolarising Noise on the 
Chain 



In [7] , the optimality of the PPT condition for detect- 
ing entanglement was proved for a vast parameter regime, 
which simply included thermal noise as a special case. 
Any instance in which all the coefficients s y were pos- 
itive was also covered. As such, we might ask how far 
beyond the thermal state our results apply, even just for 
chains. This is certainly an interesting avenue for future 
study, however proofs are going to be vastly more difficult 
than for GHZ states. For instance, in the full separabil- 
ity proof, for star graphs, if K x is compatible with K y 
and K z , then K y and K z are compatible, and this makes 
it very easy to group terms together. However, this is 
not the case for graphs such as the chain (consider, for 
instance, K%, K\Ka and K\K^). Still, we can briefly con- 
sider how likely it seems that this should hold. Do the 
PPT and full separability thresholds coincide? Are the 
optimal eigenvector and bipartition the same as for the 
thermal state? Do the NPT and distillability thresholds 
coincide? 

We shall provide some hints to these answers by con- 
sidering the graph-diagonal state 



P : 



|V> {i>\ 



which we will again address for the linear graphs. This 
noise model was considered in [3], where it was observed 
numerically, using the distillation protocol of [TBI HSL 
that distillation of the state was possible at a > 2 L N / 2 l ( 
but the impossibility of distillation was only proven for 
a < 2. We will analytically show that for even N, the 
distillation protocol of [TH] [19] does indeed distil for all 
a > 2 N / 2 , but we shall start by showing the converse, 
that distillation is impossible for a < 2 L^/ 2 ] , and our 
particular reason for doing so is to illustrate that the 
choices of x and z that are optimal for the thermal state 
are not optimal in this case. 

The eigenvalues of p, indexed by x, under bipartition 
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z are given, up to normalisation, by 

2 N + a (-l) :c ' 3 '(-l) E « Sfe?M + l( *" eZn + l) , 
ye{o,i} N 

where the sum is just /^(l). Let's consider taking the 
two-colouring bipartition of the chain, and x = {11 ... 1}. 
Using the results of Sec. |III[ we can evaluate this and find 
that the state cannot be distilled if 

2 N + 2 (JV+l)/2 a CQS ^ N + 1) / 4 ) > 0. 

Hence, for N mod 8 = 2, 3, 4, we have that the critical 
value is a — 2>- N / 2 l . For other values of N mod 8, there 
are other choices of a: which give this result. For instance, 
for N mod 8 = 4, 5, 6, the same value is achieved using 
a; = {11... 100}. 

Now let us instead consider z = {100. . .0} and x — 
{1100 ... 0}. The corresponding eigenvalue remains neg- 
ative until a < 2, so the optimal choice of z for prob- 
ing PPT across all bipartitions is certainly not the two- 
colourable bipartition. We have also verified by brute- 
force that for N < 6, the state is indeed fully separable 
at a = 2. 



1. Optimal Distillation 

Let's consider a two-colourable graph G, with the two- 
colourable bipartition dividing the vertices into sets Va 



and Vb of sizes Na and Nb respectively, N = Na + Nb ■ 
We are interested in the distillation of a state 

_ t + a\^) (ip\ 
P ~ 2 N + a ' 

which has diagonal elements X x = (A^plA^) = (1 + 
<x$x,o) I (2 W + a). For convenience, we will split the string 
x into fJ-AfJ-B with fj, A G {0, 1} Na , \i b £ {0,1}^. The 
distillation protocol of (THUS] gives two protocols PI and 
P2 which arc defined in terms of their recursion relations 
for these diagonal elements, 

pi . \ (n+l) = A (n) A (") 

P2 ■ A (n+1) =Vl("> A (n) 

If all the coefficients initially obey 

Ao^b = A to V/iB 7^ 

A MA = A-rO V/iA 7^ 
AjUjtjUB = A^ VfJ.A,VB 7^ 0; 

as they do for the globally depolarised state, then after 
the application of PI or P2, the new coefficients also 
have this same structure. One can write down the action 
of PI followed by P2. 



A^ +2) = (A^ 2 + (2^-l)A£ )2 ) 2 + (2^-l)(Ai" )2 + (2^-l)Al«) 2 ) 2 

A£ +2) = A^ )2 (2A&> + (2^ - 2)A^) 2 + (2 N - l)X^f (n® + (2 N * - 2)A£>) 2 

A£ +2) = (a^ )2 + (2^ ~ l)Ai"/) (u®* + 2(2"- - 1)A^ 2 + (2^ - 2)(X^f + (2^ - 1)A£"/)) 

A!"/ 2 ) = (4A^ + 2(2 W - - 2)A&> + (2 N * - 2) ( 2 A$ + ^ B ~ 2)A&>)) (2A$ + (2 W * - 2)A^) 



Under the assumption that Aoo _ A^Oj Aox, ^xx _ 0, and 
recalling that the coefficients must be normalised such 
that 

A 00 + (2 Na ~ l)A a0 + (2 Nb - l)Ate 

+ (2 W - 2 ^ - 2^ + 1)A M = 1, 

we find that there are two trivial fixed points of the map, 
a A A A 1 — / {MAO} 

1^00, AaO, Aoxj A^j- — < {111 1}/2 JV 

These are stable fixed points i.e. local attractors. 

If we select N even with Na — N/2, as is the case 
with a linear graph, then a third trivial point can be 



found with Aoo = 1/2 Na for all possible values of A^o, Aoa; 
and X xx . This is an unstable fixed point. These are the 
only valid fixed points subject to the aforementioned con- 
straints, and hence we conclude that repeated alternate 
applications of PI and P2 serve to distil the state pro- 
vided the initial state has A o > 1/2 Na , i.e. 

1 + a 1 

2 N + a > 2^2 

or a > 2 N / 2 , analytically proving the numerical obser- 
vations made in [3], but also proving optimality of this 
solution. 
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Appendix B: Local Depolarising Noise on the Chain 

Even if the choices of x and z which are optimal for 
the thermal state are not optimal for the fully depolaris- 
ing case, the existence of a PPT bipartition is necessary 
and sufficient for the impossibility of distillation, and, 
certainly for small sizes, N < 6, the existence of an NPT 
bipartition is necessary and sufficient for the persistence 
of entanglement. Let us now treat a different case, that 
of local depolarising noise 



3p 
4 



P 



(X n pX n + Y n pY n + Z n pZ n ) . 



For a general bipartite graph with adjacency matrix A, 
we have that 



Wy+WA y mod 2—y-(A-y mod 2) 



For the special case of a chain of N = 4 qubits, we 
can evaluate f x ,z{p) by brute force and find the optimal 
choices of x = {1100} and z = {1000}, corresponding to 
a critical point described by 

l-4(l-p) 3 -5(l-p) 4 = 0, 

yielding a maximum value of p ~ 0.468 at which the state 
becomes PPT with respect to all bipartitions. In compar- 
ison, our separable decomposition shows that the state is 



separable if p > 1 - \/(2v / 2 - l)/7 « 0.489. It appears 
that the two do not coincide, and so we should not expect 
the PPT condition to be as universal for the linear graph 
state as it was for the GHZ states. The techniques of [M] 
failed to witness any of this entanglement, but this was 
limited by the finite computational resources available, 
and is hence inconclusive. 



Appendix C: Thermal Graphs 

Ideally, we would like to make similar studies to those 
of the tree graphs for general thermal graphs, i.e. to be 
able to calculate the critical temperature up to which the 
PPT condition is violated, and to show that beyond that 
threshold, the state is fully separable. To date, the com- 
binatoric factors have prevented us from giving a proof, 
but for completeness, we present those considerations 
which are possible. 

Although we do not know the optimal bipartition 
for non-two-colourable graphs, evaluation of the partial 
transpose condition of G will always be identical to the 
evaluation on a bipartite subgraph of G. Hence, it suf- 
fices to restrict to considering two-colourable graphs. 

We have already seen how one can evaluate the critical 
temperature for GHZ states [7] and linear cluster states. 
We would like to do the same for a 2D graph in spite of 
any unproven hints from computational complexity. To 
date, we have been unsuccessful, and have been forced to 



u 3 



M, for an MxM square lattice 

4 5 fT 



Distillability Bound 

FIG. 3: Plot of the critical temperature for an N = M x M 
square lattice. Each point is a lower bound on the critical 
temperature for larger M. 

resort to numerics. The difficulty here is that to evalu- 
ate f(s), we need to sum over an exponential number of 
terms. Some advantage can be gained by rewriting the 
adjacency matrix of a two-colourable graph in a block 
structure 



1 ' A tc 



where A tc is an Na x Nb matrix, and we will assume 
that N A < N B . This yields 



E 

xe{o,i} r 



\N B — w al .A t 



. ( C1 ) 

providing a sufficient speed-up to enable some basic nu- 
merics on a square lattice of N — M 2 qubits. We have 
plotted the critical temperature in Fig. |3j which suggest 
a tendency to a finite critical temperature for large N, 
but certainly does not prove it. 

A further efficiency saving can be made by observing 
that the set of strings x that satisfy x ■ A tc mod 2 = 
form a subgroup of the group of binary strings under 
addition mod ulo 2. It is only necessary to perform the 
sum in Eqn. (CI I over the strings x £ {0, 1} Na that are 



in the coset to the subgroup. 



1. Full Separability 

The other extreme to showing the presence of entangle- 
ment is its absence, i.e. the full separability of the state. 
We have already seen how for trees this corresponds ex- 
actly to the PPT threshold. We report our progress to- 
wards a similar proof for two-colourable graphs. 

such that 



Let us define h x = (-l)™*(-l)£{i,i}eE 



/(*: 



If y is compatible with x, then h x = 



hyh x( $ y because there must be an even number of edges 
between the two blocks (we will prove this in Lemma |2|. 

We will now specify how to make our decomposition. 
For each K x , find a compatible decomposition (and de- 
compose each term as far as possible). Let's assume that 
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K x decomposes into K y and K x q v . Let's also further which can be rewritten as 



assume that the coefficients s y and s x ^ y appear in the 
1 term as h y s y and h x ^ y s xl $ y . We want to see that s x 
appears as h x s x . We use the decomposition 



s x (t - h y K y )(t - h y K x(By ) 

This means that s y i-> s y + h y s x and 
so that, in the 1 term, they become 



I ~\- hyS X 



hySy 



1 + hyh X ($y) 



Since h y = l and h y h X (£, y — h x , this certainly gives h y s y + 
h X (8 y s x( $ y +h x s XJ as desired. This inductively proves that 
all the terms in front of the 1 term are exactly those 
that appear in /(s), assuming that the first term of the 
induction can be proved i.e. every term K x that does not 
have a compatible decomposition appears as h x s x in the 
coefficient in front of the 1 term. However, we know that 
if there is no compatible decomposition, we are using 
s x (t + K x ), and therefore, s x appears in f(s) as — s x . 
Thus, if our decomposition is to hold, we require that 

Lemma 4. If h x — 1 then K x has a compatible decom- 
position. 

Proof. We start by translating the "compatible basis" 
decomposition to graph theoretic terms. Consider x 
as a bit string with bits corresponding to the vertices 
of the graph, y is a compatible decomposition pro- 
vided yi — if Xi — and, for all sites i, either (1) 
Xi = and J2{ij}eE x j mod 2 = or (2) y t = 
and E{ij}e_B% mod 2 = or (3) (x - y)i = and 
Ylu j}£E^ x ~ v)o mod 2 = 0. Note that if Xi = 0, then 
the only condition that we have is yi = - it does not 
matter what the neighbours are doing. For that reason, 
it suffices to consider only induced subgraphs G' of G 
defined by Xi = 1 (note that if G' is not connected, there 
is a trivial decomposition given by the distinct graphs). 
Hence, since h x — 1, G" has the property that the num- 
ber of edges + number of vertices is even, and y and x®y 
represent a non-trivial partitioning of the vertices V of 
G' into Vi and V 2 , Vi U V 2 = V . 
At each site i, we either have 

yi = and yj mod 2 = 

or 

(x — y)i = and (x — y)j mod 2 = 0, 

where A is the adjacency matrix (edge set) of G' . These 
two terms are mutually exclusive, and we can hence write 
penalty terms to evaluate when these are violated, 

Y^(l- yi )(l-(Ay mod 2),) 

i 

+Vi(l - {A(x -y) mod 2)i) = x T x, 



(x - y) T (Ay mod 2) + y T (A(x - y) mod 2) = 0. 

Since both terms are non-negative, this means that both 
must be 0. The interpretation of this equation is that 
each vertex in either V\ or V% must have an even number 
of edges crossing the partition. Thus, the corresponding 
graph theoretic lemma is 

Lemma 5. For a connected two- colourable graph G' with 
(number of edges + number of vertices) even, there are 
two non-empty vertex sets V\ and 2, such that each vertex 
in V\ has an even number of neighbours in V2 and vice 
versa. 

Lemma [5] is exactly the property that is referred to as 
the existence of an 'Eulerian Edge Cut', and its existence 
has indeed been proven for a two-colourable graph when 
(number of edges + number of vertices) is even [35] 3 • □ 

This proves that the coefficient in front of the 1 term 
in the separable decomposition is /(s)/2 for a two- 
colourable graph, and is hence zero at the critical PPT 
temperature, becoming negative above that temperature 
(i.e. not a valid decomposition). This is true for all 
two-colourable graph states (and thus, graph states of 
all graphs which are equivalent to two-colourable ones 
under local unitary rotations). Now, however, the chal- 
lenge is to prove that the coefficients in front of the non-1 
terms are also positive, which would prove that beyond 
the PPT threshold, a thermal graph state is fully sepa- 
rable. To date, we have only achieved this feat for the 
special cases of the GHZ state and the thermal tree graph 
states. We have also verified this to be true for all even 
rings of N < 120 qubits (the coefficients are readily re- 
lated to those of a chain, for which we gave the solution 
in the main text). Indeed, if evaluating /(s) for a given s 
turns out to be a hard computational problem (Sec. II A I , 
one would also expect the same to be true for the differ- 
ence between f(s) and the coefficients in front of the K x 
terms. 



Appendix D: Separability of Star Graphs 

In [7], we gave an identical form of separable decom- 
position to that presented here, except specialised to the 
case of star graphs. This had the massive advantage that 
all compatible terms are mutually compatible, and this 
enabled us to give a very broad condition on when the 
PPT threshold coincided with full separability. Never- 
theless, there were still examples of GHZ-diagonal states 
for which the two did not coincide. It was observed that, 
at least for special cases, it can be proven that there is 
entanglement beyond the PPT threshold, but also that 



We thank Simone Severini for pointing out this reference. 
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the constructed separable decomposition was not opti- 
mal. This last step was done by magically coming up 
with a new separable decomposition. Here we want to 
illustrate the logical formulation of how we did this so 
that it can be more widely applied. 

Our starting point is a star graph, where vertex 1 is the 
central node (root) connected to N — 1 leaves. We know 
that for x, z e {0, l}^ -1 , the PPT condition is tested by 
positivity of all the 

f X0Xfi M = E {-i) x - y s ay +{-ir{-i) {xmzyv siy 

yefo.i}"- 1 

On the other hand, we can give a fully separable decom- 
position, based on the previously expounded technique 
of finding a compatible basis, of the form 

P = E l*ivl(l + sgn(*i,,)#i,,) 
j/efo.i}"-! 

+ (j2 s »v K »y 1 xe ^ N -, E *°«(-ir v ) 

\ y v / 

+ 1 L^E^Mr^-EKi)- 

V v y I 

This state is positive provided 



mm 



E 

yefo,!}"- 1 



■S0y(-iy 




> 0, 



and hence this is easily matched with f Xo x,0z(s) over a 
large range. For N = 3, the condition that needs to be 
satisfied is 



n s i^° 

y£{0A} 2 



(Dl) 



A state such as 



Now we will concentrate on the K\ y terms, and try to 
give a better decomposition of them. The first observa- 
tion that we need to make is that if we had two strings 
of tensor products of Pauli operators [pc 7^ fu) 

SABC&A ®<JB®°C + SABD&A ® <7 B <8> <7l5, 

then in our previous decomposition, we would have added 
(sabc+ s abd)^ to make these positive. However, we can 
rewrite them as 

<ja®p b ® {sabc°c + sabd^d), 



which only requires \/s\ BC + s 2 ABD t to make it posi- 
tive. In fact, for a graph state, this situation can never 
arise since any two products of stabilizers always differ 
on at least two sites. However, we will now see how to 
'twist' the Pauli basis such that it can potentially arise. 
Note, however, that this technique is always going to in- 
troduce square roots, which are never present in f(s), so 
this technique, while having the potential to improve on 
the previous separable decomposition in some parameter 
regimes, can never widen the regime where PPT and full 
separability coincide. The next step is to observe that for 
a graph where there is a vertex of degree 1, then there are 
stabilizer terms composed of only 2 Paulis, which means 
that products differ by only two Paulis. For instance, in 
the three qubit case we have 



x®z®z 

Y®Y®Z 



-Y®Z®Y 
-X®Y®Y 



(note that one has to be very careful with any -ve signs 
that may appear from multiplying out the stabilizers). 
We can expand the X and Y on qubit 1 using \{X + 
Y) ± ^(X — Y). The fact that there are products of 
stabilizers which differ by only two Paulis, one of which 
is the spin we've just expanded, means that there will 
now be terms that differ on only one site each, 



-\{X + Y) 



8(1 



/ 3 \ 5 \A +1)1 

— ) \T[(l + K n )-2K 1 K 3 + al\. -\{X + Y) 



Z^ 
)Z< 

i"(X + Y) ® Y $ 



z 

i)Y 
Z 

OY 



l(X-Y), 



Z( 
Z( 



\{X-Y)®Y 
-I(X-Y) 



Y 



i Z 
Y 

%>Z 
R) Y 



does not satisfy this condition. For a > 2, p is a valid 
state, but also PPT. On the other hand, the above sepa- 
rable decomposition only functions if a > 4. Our aim is 
now to show constructively how to improve this bound 
up to the threshold of a = 2^/2. The first thing to ob- 
serve is that the grouping of all the compatible terms is 
probably as efficient as it's going to get so, as before, we 
will use 

(e s °v K °y 1 xe{ ™if»-. E M-ir^ 

which, in this case, is just (1 + K 2 )(l + K 3 ). 



and we can collect them as above (we could either collect 
on qubit 2 or qubit 3), 

\{X + Y)®Z® {Z-Y) \{X-Y)®Z®{Z + Y) 
\{X + Y)®Y® (Z-Y) -\{X-Y)®Y®{Z + Y)- 

Now we see that these terms only differ on one site (this 
is because there is a vertex which has two neighbours 
with degree 1), and so can also be combined 

\{X+Y)®(Z+Y)®{Z-Y) + \(X-Y)®(Z-Y)®(Z+Y). 

Each of these terms has eigenvalues ±\^2 and hence 
a = 2\/2 suffices to give a separable decomposition (we 
also succeeded in witnessing the entanglement between 
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2 < a < 2\/2 using the numerical techniques of [33]. Us- 
ing a modification of [5J, it can be proven that the state 
is entangled up to a — 2\[2 [36 . This rearrangement of 
terms to form a new separable decomposition will always 
give some advantage whenever Eqn. (Dl) is not satisfied 



for the 3-qubit star graph (although the general expres- 
sion is very messy to write down), and can presumably 
be applied in other graphs with star-like protrusions if 
the signs of the weights of s y happen to be of the correct 
form. 
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